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We consider the possibility that classical dynamical systems display motion in their lowest energy 
state, forming a time analogue of crystalline spatial order. Challenges facing that idea are identified 
and overcome. We display arbitrary orbits of an angular variable as lowest-energy trajectories for 
nonsingular Lagrangian systems. Dynamics within orbits of broken symmetry provide a natural 
arena for formation of time crystals. We exhibit models of that kind, including a model with 
traveling density waves. 



In this paper we will investigate a cluster of issues 
around the question of whether time- independent, con- 
servative classical or semiclassical systems might exhibit 
motion in their lowest energy states. Fully quantum sys- 
tems are the subject of a companion paper [I]. Related 
issues have been raised in a cosmological context 
but those investigations consider quite diflFerent aspects, 
in which the time dependence introduced by the expan- 
sion of the universe plays a significant role. 

General considerations. When a physical solution of a 
set of equations displays less symmetry than the equa- 
tions themselves, we say the symmetry is spontaneously 
broken by that solution. Here the meaning of "physical 
solution" can be interpreted differently in different con- 
texts, but one interesting case, that will concern us here, 
is of the lowest energy solutions of a time- independent, 
conservative, classical dynamical system. If such a solu- 
tion exhibits motion, we will have broken time transla- 
tion symmetry spontaneously. If the dynamical variable 
is an angular variable, then the motion will be periodic 
in time, so the time-translation symmetry is not entirely 
lost, but only reduced to a discrete subgroup. Spatial 
periodicity is, of course, associated with formation of or- 
dinary crystals, so it is natural and suggestive to refer to 
the formation of time crystals. 

It is very easy to construct simple Lagrangians or 
Hamiltonians whose lowest energy state is a spatial crys- 
tal. With 4>{x) an angular variable, the potential energy 
functions 
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with all the Greek coefhcients positive, are minimized for 
^ = fl' ^ = ^\f^ respectively. In both cases the 
spatial translation symmetry of the original potential is 
spontaneously broken; in the second case inversion sym- 
metry is broken as well. The combined inversion 4'{x) — >■ 



— 0(— a;) is preserved in both cases, as is a combined in- 
ternal space-real space translation (j)(x) cj){x + e) — ^e. 

From this one might surmise that time crystals are like- 
wise easy to construct, at least mathematically. On sec- 
ond thought, however, reasons for doubt appear. Speak- 
ing broadly, what we're looking for seems perilously close 
to perpetual motion. Also, if the dynamical equations 
conserve energy, then a minimum-energy solution where 
the variables trace out an orbit implies that the energy 
function assumes its minimum value on a whole curve in 
{(j), 4>) space - not, as we expect generically, at an isolated 
point. 

Dynamical equations. That easy/impossible di- 
chotomy carries over into the dynamical equations. If 
one simply turns the space derivatives in Eqn. ([I]) into 
time derivatives, then the resulting Lagrangians 
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are associated with the energy functions 
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0, the second at 



jj^. So the analogue of our first symmetry- 



The first of these is minimized at 

breaking example in Eqn.Q has collapsed, but the sec- 
ond survives, with a quantitative change. 

On the other hand if we convert the space derivatives 
in Eqn. ([T]) into momenta, the resulting Hamiltonians 
are 
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We find precisely the original algebraic structure for the 
minimum-energy solutions, viz. pi — j^,P2 = -'-^/ff 
respectively. Their physical implications are entirely dif- 
ferent, though. Indeed, they correspond to 0i = 02 = 0: 
thus no symmetry breaking occurs, in either case. 

This disappointing consequence of the Hamiltonian 
formalism is quite general. The Hamiltonian equations 
of motion 
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indicate that the energy function i?(pj(0), g-' (O)) = 
H{pj{Q), (0)), regarded as a function of the dynami- 
cal variables at a chosen initial time, is minimized for 
trajectory with pj — — 0, since the gradients on the 
right-hand side of Eqn. ^ vanish. 

How do we reconcile this very general null result in 
the Hamiltonian approach, with our positive result in 
the Lagrangian approach? The point is that the La- 
grangian L2, which gave symmetry breaking, cannot be 
converted into a Hamiltonian smoothly. Indeed (putting 
K2 = K, A2 = 1 for simplicity) expressing the algebraic 
recipe for the Hamiltonian 
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as a function of 
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leads to a multi-valued function, with cusps where ^ 



T^fl72~j corresponding precisely to the energy 



0, i.e. p ~ 

minima p = ^V^' (^^^ Figure 1.) For k < H{p) is 
regular, but as n passes through zero there is a swallow- 
tail catastrophe. 

At the cusps the usual condition that the gradient 
should vanish at a minimum does not apply, and so our 
null result for smooth Hamiltonian systems is avoided. 

For classical physics the Lagrangian formalism is ad- 
equate, so let us follow that direction out further. A 
logical next step would be to add a potential 1^(0) to i2- 
Doing that, however, leads us directly into the problem 
with energy conservation we anticipated earlier. Mini- 
mizing V ^ we will find a preferred value for (p — (pQ. But 
minimizing the kinetic part will favor motion in 0, and 
there is a conflict. 

We can elucidate this issue as it arises for a general 
Lagrangian system. Suppose that the energy function of 
a system with many degrees of freedom is minimized by 



FIG. 1: Energy is a multivalued function of momentum. 



nonzero velocities (pQ 7^ 0, so that 
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Then in the equations of motion 
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the coefficient of the acceleration in the direction (p^ oc 0g 
vanishes at 0q. This leads to a difficulty if there are 
forces that require acceleration in that direction. The 
equations of motion, which generally serve to determine 
the accelerations, may not admit a consistent solution if 
the coefficient of the accelerations cannot be inverted. 

As we shall discuss below, there are physically inter- 
esting models that avoid any singularity. Before turning 
to those, however, we now discuss a prescription for con- 
tinuing through it. 

Brick Wall Solutions: Upon integrating 



directly we obtain 
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where the ± signs are independent. 

The argument of the inner square root is non-negative 
if and only if 
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where A = E Eq > is the energy above the minimum 
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kinetic energy Eq = — fj- The inequahty is saturated 
when (j) = zL^y~^, i.e., when the kinetic energy is mini- 
mized. Close to a point cjjt where this happens. 
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Since cannot continue past 0t without violating the 
bound (12), it suddenly reverses direction, </> = ^\/~^ ^ 



Such reversal conserves energy, but requires a 
sudden jump in the momentum. A "brick-waU" turn- 
ing point is appropriate to a potential enforced by an in- 
finitely massive source. Unless 0t is an extremum of V((f>) 
the acceleration diverges at (j)t^ as anticipated above. 

Small oscillations about the minimum of a generic po- 
tential V{(j)) ~ ^m(0 ~ 0o)^ exhibit turning points of 
this type, with bounded orbits that oscillate between 
0t = 00 — v^SA/Tt and 4>o + y^2 A/ fi . In the limit of 
small A, the orbits ricochet about the minimum, with 
nearly constant speed \(f)\ — -\/^, reconciling the appar- 



ently contradictory conditions (j) = i^y^ and 4> = 4>o. 

More con ventional turning points arise if the inner ± 
sign in (111 is negative, i.e. for V{4)) > E. Now (pt ~ 
at turning points where V{(j)t) ~ E, and the particle 
changes direction smoothly. 

fgh model: Lagrangians of the form 
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for functions f {(()), g{(j)), h{(f)) lead to energies of the form 
E = 3f^'^ + g^'^-h 
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(Note that / may be absorbed into a redefinition of ip; 
then for constant g this reduces to the model of Eqn. 



(101.) If 
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the energy will be minimized along the curve (fP' — 
0, for any / > 0, g < 0. Thus we have solutions 
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This construction demonstrates that any orbit with a ve- 
locity that does not change sign can be realized, in many 
ways, as the stable minimum energy solution to an ap- 
propriate, reasonably simple Lagrangian. It would be 
interesting to explore the (thought-)engineering of more 
complex models of this kind, and in particular, to in- 
vestigate the possibility of realizing non-dissipative ( "re- 
versible") computation in the lowest energy state of a 
dynamical system [1]. 



Choosing the constant in ( 16 1 to be zero, constancy of 
the energy E >Q leads to 
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This equation is of a familiar form; it expresses the con- 
servation of a pseudo-energy E for a particle with mass 
m = I and a two-branched i?-dependent potential, ac- 
cording to 
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This result allows us to infer the qualitative dynamics, 
based on familiar mechanical concepts. Perhaps the most 
interesting question is the existence, or not, of turning 
points. Putting = into Eqn. ( 19 ) we find that 
E — —g^{(f>t)/Gf{4't) = Knax evaluated at the turning 
point (s) (pf The motion is confined to a region where 
V < l^nax- Thus the model can support motions in 
which the velocity changes sign, but these motions re- 
quire higher energy that the minimal orbit, which is 
unidirectional. Actually nothing in our analysis of this 
model has depended on the periodicity of (p; upon drop- 
ping that assumption, we find the curious situation that 
some unbounded motions have smaller energy than any 
bounded motion. 

Avoiding Singularity Semiclassically: If we relax the 
condition Eqn. (16 1, by allowing a non-constant W = 
— ^ h, wc find that there are initial conditions for 
which the equations of motion eventually become singu- 
lar, as discussed previously. We are guaranteed to avoid 



such conditions if in Eqn. (151 we require 
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(Note that in any case E > Wmin-) Thus we have models 
that work smoothly for high energy, but become singular 
at low energy. That is the opposite of the usual philoso- 
phy of effective field theory; however it does correspond 
to the use of perturbative QCD. 

We might also expect that the quantum-mechanical 
versions of these models might be more robust, in that 
the uncertainty in position and velocity might smooth 
over a small region of singularity. Genuine quantiza- 
tion of such models - or, for that matter, of the tuned 
fgh models and the natural, locked models to come - 
presents interesting issues. The Lagrangian formulation 
is adequate for path integral quantization, and the higher 
derivative terms tend to damp the contribution of the 
most irregular paths, particularly if we continue the time 
to have a small negative imaginary part. So there are 
no obvious show-stoppers, but no existence proofs ei- 
ther. The Hamiltonian formulation poses different is- 
sues. As we've seen, in interesting cases the Hamiltonian 
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is a multi-valued function of the momentum. This im- 
plies that the momentum does not provide a complete 
set of commuting observables. Nor, therefore, does the 
position. Wave functions must be defined as functions 
of expanded spaces. We will report work on this subject 
elsewhere. 

What we can discuss simply is semiclassical quan- 
tization. Thus we consider orbits obeying a Bohr- 
Sommerfeld condition 
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with n an integer, and S a correction for turning points. 
(For simplicity we specialize here to/=|,g = — |, and 
h = 0.) If we ignore, at first, the potential, then the 
minimal energy orbit is at = ±-^/^, and for it 
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If this expression is not equal to 2'!Th{n + 6), the quantiza- 
tion condition will lead us to a nearby higher-energy orbit 
for the ground state, with some \n\ 1 in the relevant 
(semiclassical) limit. If the potential is small enough, 
that extra energy will be enough to enforce Eqn. (21), 



and keep us out of the region where the equation of mo- 
tion breaks down. Wave packets constructed from n near 
the preferred value will describe, approximately, the mo- 
tion prescribed by the classical dynamics. 

Naturally Flat Directions; Double Sombrero: It can be 
natural to have energy constant along an orbit, if the 
points of the orbit are related by symmetry. If we want 
this situation to occur along a trajectory representing 
non-trivial motion in the minimum-energy state, then 
the points assumed at different times must be related 
by symmetry transformations, which implies that none 
of them is invariant. So we will be looking at models 
with spontaneously broken symmetry. 

Consider first a Lagrangian with a "sombrero" kinetic 
term: 
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The matrix of second derivatives of L with respect to tp^^ 
appearing in Eqn. Q is 
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This has a zero eigenvalue with eigenvector 

cisely when = + -01 = k/3. 

If the potential V has a one-parameter family of degen- 
erate minima, the minimum-energy solution will move 
along the trough of V at constant speed k/3. The 
potential 
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is symmetric under tjji-'ip2 rotations and, combined with 
the kinetic term in Eqn. (23), defines a "double som- 



brero" model, with circular motion at constant speed in 
the lowest-energy state. 

Alternatively we may rewrite this model and its gen- 
eralizations in terms of polar fields p and (j), where 
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if. Then the double sombrero La- 



grangian takes the form 
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If p is set equal to its value ■\/2/i/A at the minimum of 
V{p), this reduces to our original Lagrangian (H). Gen- 
eralizing, any Lagrangian with a kinetic term that is a 
polynomial in (f), p, and p, and a potential energy depend- 
ing only on p, will preserve the symmetry (f> (f> + rj. 

Charge and Locking: The charge operator associ- 
ated with the original (broken) symmetry is Q = 
— J i{ip*7ript — i-piXif) where TTtp = ^ depends only on 
and p. Thus in states with constant, non- vanishing val- 
ues of p and we have a non-zero, uniform density of Q. 
This is significant in two ways: 

First: If we suppose that our system is embedded in a 
larger symmetry-conserving bath and undergoes a tran- 
sition to the symmetry-breaking state, e.g. that it is a 
material body cooled through a phase transition, then 
the transition will necessarily be accompanied by radia- 
tion of an appropriate balancing charge. 

Second: Although invariance under both infinitesimal 
time-translation 0(t) — >■ <^(t -f- e) and under infinitesimal 
phase (charge) translation — > + 77 are broken by solu- 
tions of the type 0(i) = uit + P - constant - the com- 
bined transformation with uje + 7] = leaves the solution 
invariant. Thus there is a residual "locked" symmetry. 
To exploit it, we can go to a sort of rotating frame, by 
using 



H = H-ujQ 



(27) 



to compute the evolution. (Here we measure charge us- 
ing the charge of as the unit.) In the rotating frame, 
the equations of motion will not contain any explicit time 
dependence, but there will be a sort of effective chemical 
potential (associated however with a broken symmetry). 
The most interesting effects will arise at interfaces be- 
tween the locked phase and the normal phase, or between 
different locked phases, as exemplified in the preceding 
paragraph. 

Space-Time Structure; More Complex States: We can 
also contemplate slightly more complex examples, that 
support qualitatively different, richer physical effects. If 
there is a potential for Viy9, or ultimately for Vp, that fa- 
vors gradients, then we can have a competition between 
the energetic desirability of putting p at the energetic 
minimum and accommodating non-zero gradients. Un- 
like in the case of time derivatives, there is no general 
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barrier to reaching a stable compromise. To keep things 
simple, let us suppress the underlying Lp structure and 
consider the potential 



ap 




with a, > 0. This is minimized by 
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This po reduces the translation symmetry to a discrete 
subgroup. Constant produces a charge density wave. 
If in addition we have a term of the form 
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then at the minimum 4'q{x) will develop spatial structure 
as well, according to 4)o{x) = ppoix) + (3, breaking the 
phase (charge) symmetry completely. (Note that V^radiont 
respects the symmetry 4> ^ (j) + rj.) 

We can engineer similar phenomena involving most 
easily if we work at the level of the energy function. One 
can derive general energy functions involving powers of 4> 
from Lagrangians of the same kind, so long as there are 
no terms linear in 6. Thus if we have an additional term 
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then at the minimum we have 
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4>o{x,t) = ppo{x,t)^p (32) 
po(a;,t) — ^J- sm(^J~^{x ± vt) + (33) 



( 33 1 wc have adapted Eqn. ( 29 1 by taking 



Here in Eqn. 

a = ±vt + a. In doing this we assume that the energy 
intrinsically associated with time derivatives of p van- 
ishes (or that it is dominated by the locking effects of 
Eqns. (|30 31)). Both spatial and time translation are 
spontaneously broken, as is reflected in the disposable 
constants a, f3, and so is time-reversal T, as reflected in 
the disposable sign. 



Combining Eqns. (32 33), we now have a traveling 



charge density wave. Thus this example exhibits its time- 
dependence in a physically tangible form. The residual 
continuous symmetry is reduced to a combined discrete 
time-space-charge transformation. Although our con- 
struction has been specific and opportunistic, it serves to 
establish the existence of a universality class that, since 
it is characterized by symmetry, should be robust. It is 
noteworthy that cyclic motion of (j) in internal space has 
given rise to linear motion in physical space. 

Comment: All of our constructions above have been 
nonrelativistic. In a relativistic theory there are relations 
among the coefficients of time and space gradient terms. 
The relativistic quartic term L cx ((9o0)^ — (V(/))^)^ leads 
to an energy that is unbounded below, for large gradi- 
ents of one kind or another. But use of a sextic en- 
ables positive energy. Indeed, the energy function for 
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For n odd this is semi-positive definite, with a zero at 
(docj))'^ — iyipY unless n = 1. For n even it has no 
definite sign. Bounded energy requires that the leading 
term has odd n and a positive coefficient, and that the 
coefhcient of the n = \ term is non-negative. This con- 
sideration might help to constrain the models of [2] [3]. 
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